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Abstract: This paper deals with the problem of trajectory tracking of a class of bilinear 
systems with time-varying measurable disturbance, namely, systems of the form x{t) = [A + 
+ d{t). A set of matrices {A,Bi} has been identified, via a linear matrix 
inequality, for which it is possible to ensure global tracking of (admissible, differentiable) 
trajectories with a simple linear time-varying PI controller. Instrumental to establish the result 
is the construction of an output signal with respect to which the incremental model is passive. 
The result is applied to the boost and the modular multilevel converter for which experimental 
results are given. 
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1. INTRODUCTION 


Bilinear systems are a class of nonlinear systems that 
describe a broad variety of physical and biological phe¬ 
nomena Mohler (2003) serving, sometimes, as a natural 
simplification of more complex nonlinear systems. There 
is an amount of literature devoted to the study of the 
intrinsic properties or to stabilization of equilibrium points 
for these systems, see for example Elliot (2009). However, 
to the best of our knowledge, there is no general result 
for the design of controllers that ensure global tracking of 
(admissible, differentiable) trajectories. 

The main objective of this paper is to provide a theo¬ 
retical framework—based on the property of passivity of 
the incremental model—to establish such a result. Our 
motivation to pursue a passivity framework is twofold, 
on one hand, it encompasses a large class of physical 
systems. On the other hand, it naturally leads to the design 
of PI controllers, which are known to be simple, robust 
and widely accepted by practitioners. Our main result is 
an extension, to the problem of tracking trajectories, of 
Jayawardhana et al. (2007); Sanders and Verghese (1992) 
that treat the regulation case. See also Castanos et al. 
(2009) for its application to PI stabilization of RLC cir¬ 
cuits and Hernandez-Gomez et al. (2010) where the result 
is used in power converters. 


An important motivation for our research is the derivation 
of simple tracking controllers for power converters. In clas¬ 
sical applications of power converters the control objective 
is to regulate the output voltage (or current) around some 
constant desired value. Modern applications, on the other 
hand, are concerned with the more demanding specifica¬ 
tion of ensuring an effective transfer of power between the 
sources and the loads—an objective that translates into 
the task of tracking time-varying references. 

In the paper our theoretical result has been illustrated 
with the application to two important problems arising 
in power electronic systems. The first one is the problem 
of Power Factor Compensation (PEG), which arises in 
renewable energy and motor drive systems with stringent 
specifications on efficiency, harmonic distortion and volt¬ 
age regulation Fadili et al. (2012); Hussain et al. (2011); 
Mather and Maksimovic (2011); Cimini et al. (2013). The 
second problem is related to High Voltage Direct Current 
(HVDC) transmission, which has recently attracted a lot 
of interest Flourentzou et al. (2009); Bahrman and John¬ 
son (2007); Ahmed et al. (2011); Zonetti et al. (2014). 
HVDC transmission consists of a grid comprising mostly 
DC lines that integrates, via voltage source converters, 
renewable energies from distant locations. One of the most 
promising candidates to integrate the topology of the grid 
is the Modular Multilevel Converter (MMC) Glinka and 
Marquardt (2003), which has several advantages with re- 



spect to its predecessors, such as high modularity, scala¬ 
bility and lower losses. In view of its complicated topology 
and operating regimes, controlling the MMC is no simple 
task. In particular, to exploit the full potential of the MMC 
it seems to be necessary to develop control strategies for 
the system operating in the rotating [ahc) frame Bergna 
et al. (2013), this is contrast with the classical strategies 
developed in fixed (dgO) frames Tu et al. (2011); Bergna 
et al. (2012). This situation leads to a tracking problem 
instead of the typical regulation one. The tracking problem 
of bilinear systems has been addressed within the context 
of switched power converters. In 01m et al. (2011) a 
methodology to track periodic signals for non-minimum 
phase boost converters based on a stable inversion of the 
internal dynamics taking the normal form of an Abel ordi¬ 
nary differential equation was presented—see also Fossas 
and 01m (2009). There are also schemes involving sliding 
mode control, for example Fossas and 01m (1994); Biel 
et al. (2004), and references therein. In Meza et al. (2008) 
the well-known passivity property of Sanders and Vergh- 
ese (1992) is used to address an “approximate” tracking 
problem for an inverter connected to a photovoltaic solar 
panel. A similar framework was studied in Meza et al. 
( 2012 ). 

The remainder of this paper is organized as follows. Prob¬ 
lem formulation is presented in Section 2. Our main the¬ 
oretical result is contained in Section 3, where a linear 
matrix inequality (LMI) condition is imposed to solve the 
tracking problem—invoking passivity theory. Section 4 is 
devoted to the synthesis of a PI controller that ensures 
tracking trajectory, under some suitable detectability as¬ 
sumptions. The result is applied in Sections 5 and 6 to 
the two power electronic applications mentioned above. 
Simulations and experimental results are included in these 
two sections. Finally, conclusions in Section 7 complete the 
paper. A preliminary version of this paper was reported in 
Cisneros et al. (2015) 

2. GLOBAL TRACKING PROBLEM 
Consider the bilinear system 

m 

x{t) = Ax{t) + d{t)+ '^Ui{t)Bix{t) (1) 

i=l 

where ^ x G M", d G M" are the state and the (measurable) 
disturbance vector, respectively, u G M'", m < n, is the 
control vector, and A G M"^", Bi G are real constant 

matrices. 

We will say that a function x* : M+ —>■ M" is an admissible 
trajectory of (1), if it is differentiable, bounded and verifies 

m 

Xi, =Ax^, + d+ ''^^u*BiX^, ( 2 ) 

i=l 

for some bounded control signal u* : M+ —>■ M™. 

The global tracking problem is to hnd, if possible, a 
dynamic state-feedback controller of the form 

z = F{x,x^,u^) (3) 

u = H{x,Xi,,Ui,), (4) 

^ For brevity, in the sequel the time argument is omitted from all 
signals. 


where F : M” x M” x —>• q G and id : M” x 

M" X —>• M™, such that all signals remain bounded and 

lim [a;(t) — x*(t)] = 0, (5) 

for all initial conditions (a;(0),z(0)) G M" x and all 
admissible trajectories. 

In this paper a set of matrices {A, Bi} has been charac¬ 
terized for which it is possible to solve the global tracking 
problem with a simple linear time-varying PI controller. 
The class is identified via the following LMI. 

Assumption 1. 3P G such that 

P = P^ >0 (6) 

sym{PA) < 0 (7) 

sym{PBi) = 0, (8) 

where the operator sym : —>■ computes the 

symmetric part of the matrix, that is 

sym{PA) = + A^ P). 

To simplify the notation in the sequel the positive semidef- 
inite matrix has been dehned 

Q := —sym(PA). (9) 


3. PASSIVITY OF THE BILINEAR INCREMENTAL 
MODEL 

Instrumental to establish the main result of the paper is 
the following lemma. 

Lemma 1. Consider the system (1) verifying the LMI of 
Assumption 1 and an admissible trajectory x*. Define the 
incremental signals 

(•):= (•)-(•)*, 

and the m-dimensional output function 

y:=C{Xif)x ( 10 ) 

where the map C : M" —)■ is defined as 

xjBj 

C{x^) := : P. (11) 

A* 

The operator ui-G y is passive with storage function 

V{x) := ^x^Px. (12) 

Hence, it verifies the dissipation inequality 
V < vJy. 

Proof. Combining (1) and (2) yields 

m in 

i={A + J2 “I” ^ ^ 

Now, the time derivative of the storage function (12) along 
the trajectories of (13) is 



Vix) 


x^P 


m m 

(A + Y. 

Z=1 2=1 


m 

= —x^Qx + UiX^ PBiXi, 
2=1 
m 

< y^ UjX^PBjX^ 

2 = 1 
m 

= y^UiX^ PBjXi, 
i=l 

= y^u, 


where (8) of Assumption 1 has been used to get the second 
identity, (7) for the first inequality, (8) again for the third 
equation and (10) for the last identity. 

Remark 1. A key step for the utilization of the previous 
result is the derivation of the desired trajectories x* and 
their corresponding control signals u*, which satisfy (2). 
As shown in the examples below this may prove to be a 
very complicated task and some approximations may be 
needed to derive them. Indeed, it is shown in 01m et al. 
(2011) that even for the simple boost converter this task 
involves the search of a stable solution of an Abel ordinary 
differential equation, which is known to be highly sensitive 
to initial conditions. 


4. A PI GLOBAL TRACKING CONTROLLER 

From Lemma 1 the next proposition follows immediately. 

Proposition 1. Consider the system (1) verifying Assump¬ 
tion 1 and an admissible trajectory x* in closed loop with 
the PI controller 


W = —x^Qx + y^u — Kiy 
= -x^Qx - y^Kpy 
< -AminlATpIlyp - < 0 

From here we conclude that the system state z,x G Coo 
and y,Q^x G £ 2 , consequently ya G £ 2 - To conclude that 
Vait) —>■ 0 it suffices, invoking Tao (1997), to prove that 
VaG C 00 . Towards this end, we first notice that x, x^ € Coo 
implies x G Coo and, this in its turn, implies from (10) 
y G Coo- Now, y,z,Ui, G Coo implies, from (14), u G Coo- 
That implies, from (13), x G Coo- Now, compute 


-1 

H- 

_ 1 


1- 

|— -K 

1_ 


Px G 


B^ 


x~^ bP 


which is bounded because x* G Coo- 


The proof of global state tracking follows noting that 
ya{t) -G 0 ensures (5) if the rank condition (16) holds. 

Remark 2. Notice that the matrix C depends on the ref¬ 
erence trajectory. Therefore, the rank condition (16) iden¬ 
tifies a class of trajectories for which global tracking is 
ensured. 

Remark 3. Condition (16) is sufficient, but not necessary 
for state convergence. Indeed, as shown in van der Schaft 
(2000), global tracking is guaranteed if ya is a detectable 
output for the closed-loop system. That is, if the following 
implication holds 

ya{t) = 0 lim [x(t) - x*(t)] = 0. 


5. APPLICATION TO POWER FACTOR 
COMPENSATION IN A BOOST CONVERTER 


z = -y 

u = - Kpy -G KiZ + u* (14) 

with output (10), (11) and Kp = > 0, Ki = Kj > 0. 

For all initial conditions (x(0),z(0)) G K" x K™ the 
trajectories of the closed-loop system are bounded and 

lim ya{t) = 0, (15) 


where the augmented output ya '- K-i- 
as 

rc(x*)i, 


im+u jg defined 


ya ■- = 


Q-^ 


with Q-^ the square root of Q given in (9). Moreover, if 

C(x*) 

. . 

then state global tracking is achieved, i.e., (5) holds. 


rank 


= n, 


(16) 


Proof. Notice that the PI controller (14) is equivalent to 
u = - Kpy -G KiZ 
z = -y. 

Propose the following radially unbounded Lyapunov func¬ 
tion candidate 

W(x, z) := V(x) + ^z^KiZ, 


5.1 Model and Problem Formulation 

The following scenario corresponds to the PEC of an AC- 
DC boost converter. Assuming linear loads and sinusoidal 
steady-state regime it is well-known that the power factor 
is optimized when the line current is in phase with line 
voltage. ^ Hence, the problem can be recast as tracking 
problem that fits the theoretical framework given in the 
previous sections. An additional requirement imposed to 
the PEC is to reduce the distortion in the current, which 
can also distort the line voltage. A standard measure 
to assess this property is the Total Harmonic Distortion 
(THD) index. 

The interleaved AC-DC boost is one of the most popular 
PEC topologies in real applications. The idea of this 
topology is to incorporate N branches with their own 
switch, each of whom admitting the control signal shifted 
by When assuming parasitic parameters such as the 
voltage drop or losses in the diodes, a typical form of 
a two branches converter is depicted in Fig. 1. Since 
these parameters are usually small, they will be neglected. 
Hence, from Fig. 1, Rjjs, Rd, Rl = 0 . Under this 
consideration, the model equations are 


^ See Garci'a-Canseco et al. (2007) for the case of nonlinear loads 
and non-sinusoidal regime. 


whose time derivative is 




Fig. 1. Boost converter circuit 

1 / /N E 

*1 = - ^(1 - M )^'C + ^ 

J^l Li 

1 E 

*2 = - ^(1 -^J-)vc + -r 
1^2 

=^( 1 + ^(1 - Ai ')*2 - ( 18 ) 

where E is the rectified AC voltage Li, L 2 ,C, R are the 
boost inductance, capacitance and load resistance, respec¬ 
tively, and fj, = 1 — /r' is the duty cycle of the switch . If it is 
assumed that Li = L 2 = ^, it follows that ii = i 2 = 
and, according to Kirchhoff’s law, = ii + i 2 (see Fig. 
1). Then, using the latter and defining x := [II and 

M = 1 — equations (18) can be expressed as 


2 2 

ii = - —ua;2 -I- yA 

Lj 1j 

11 

a;2 - ;^^2, (19) 

Clearly, the system (19) can be written in the form (1) 
with the following definitions 


, 0 0 
^ ~ 0 _ E 

RC 


fn —^ 

d= L R = V ^ 

, U „ , JD In 


A matrix P that verifies the conditions of Assumption 1 is 


V . 

Furthermore, the passive output y is 
y=xj B^x 

=Xii,X2 - X2i,Xi. ( 22 ) 

Now, convergence of the state is tested by means of the 
matrix 

■^/ \~\ ^2'k ^1* 

=00 

Q-^ 0 ^ ’ 

L vi?J 

which satisfies the rank condition (16), provided X2 is 
bounded away from zero, which is consistent with the 
operating mode of the boost, i.e., X 2 > 0. 

Reference signals a; 2 *, xi* and u* are derived below. 


5.2 Generation of references 

It is assumed that the dynamics of the output voltage X 2 
is slower respect to that of the current xi. The objective 
of control is to drive x to x* in such a way that 

(1) xi* (X E (Power Factor Correction) 

(2) Xii, (X 

(3) X 2 i, is a desired output voltage. 

Then, the resulting reference for the inductor current is 
shown in the following equation: 


" WrmsT 

where (j) is the output of a linear compensator Gcvis) 
whose input is the voltage error and Erms is the RMS 
value of the rectified AC voltage. The design of the 
compensator Gcv{s) is based on the perturbation and 
linearization technique applied in the loss-free resistor 
(LFR) boost model Erickson and Maksimovic (2001). This 
results in equivalent small-signal circuit for the design of 
the following transfer function Ridley (1989) 

^ 2 ( 5 ) _ Pgy j', 

0(s) “sCA2$’ ^ ^ 

where X 2 {s) and A 2 are respectively the Laplace small 
signal of X 2 and its DC output voltage, also (j){s) and 
$ are respectively the Laplace small signal of </> and 
its DC compensator output control signal, Pav is the 
average rectifier power and G is the boost capacitance. 
This transfer function neglects the complicating effects of 
high-frequency switching ripple, and is valid for control 
variations at frequencies sufficiently lower than the AC 
fine frequency. Ultimately, the linear compensator Gcv(s) 
consists of a PI with anti-windup technique Erickson and 
Maksimovic (2001). The output voltage controller must 
have sufficiently small gain at frequency and minimize 
the negative effect of the Right Half Plane (RHP) zero. 
Hence its bandwidth must be low. As a rule of thumb, 
setting the overall control loop bandwidth to a third 
of the RHP zero is enough to provide the closed loop 
stability. It requires, however, a compromise in the control 
performance. Erickson and Maksimovic (2001). 

The design of the controller is concluded with the calcula¬ 
tion of M*. Thus, from the first equation of (19) 

2E - Lxii, 






Table 1. 


15 


15 


Indexes 

PF 

THD 

Proportional 

96.7% 

23.2% 

Hyperbolic Tangent 

98.2% 

21.1% 


5.3 The Controller 

In summary, the resulting controller equations for the PFC 
AC-DC Boost Converter are the following. 




Simulation Trend Experimental Trend. 


Z = -y 

u = -KpU + K^z + M* 


y = Xii,X2 - X2*Xi 

2E — Lxii, 


a;i* = 


E(j) 


{Erms? 

where (f> is the output of the anti-windup PI and X 2 * is the 
desired DC output voltage. 


Remark 4- It is customary to implement a PI controller 
with changing gains. Usually nonlinear functions as the 
tanh(-) are introduced. In this case, a modification in the 
second equation of (26) leads to 

u = —6tanh -|- KiZ + m* (27) 

where a and b are free parameters to be adjusted. A 
typical shape of function 6tanh(-) is depicted in Fig. 2. 
As a matter of fact this change does not compromise the 
stability of the system. It comes down to the fact that 
the function tanh( 2 /) can be considered as the product 
kp{y)y where kp{y) is scalar function taking nonnegative 
bounded values except at y = 0 where it is zero. A practical 
comparison between controller (26) and its modification 
will be shown in the sequel. 


5-4 Simulations and Experimental Results 


Fig. 3. Inductor Current iR and Input Voltage Vg Trend. 



Fig. 4. Output Voltage vc Transient. Comparison between 
(26) and (27) 



Tinie{us) ^ . jq '^ 

Simulation Trend. 


Time{us) 

Experimental Trend. 


Fig. 5. Output Voltage vq Trend with a load change of 
30%. Comparison between (26) and (27) 


The interleaved Boost converter was simulated under 
the following considerations: the output voltage reference 
Vref = 15V, the PWM frequency PWMpreq = lOOKHz 
and the input voltage Vg = 9\sin(2Trfiinet)\ where fime = 
50Hz. 

Boost inductance, capacitance and load values are respec¬ 
tively: L = 56y,H, C = 3047/rF, R = 220. The gain for 
the compensator Gcv (s) and for the passivity PI have been 
respectively set to kp(comp) = 0.011, kn^comp) = 0.03 and 
Kp = 0.013, Ki = 0.0001 whereas the hyperbolic tangent 
parameters are a = 55 and b — 0.25. 

Figure 3 shows the trend of the inductor current ip and of 
the input voltage Vg in an interval time during the steady 
state: it is clear that the current ip (the blue line) and 
input voltage Vg (the red line) have the same shape, except 
in the instants in which the value of Vg is low driving 
the boost to the Discontinuous Conduction Mode (DCM). 
Also performance indexes, displayed in Table 1, confirm 
that the control reaches the main objective and suggest 
the use of the hyperbolic tangent instead of PI because of 
its slightly better performances. Figure 4 shows the trends 
of the output voltage vc- Its reference value, namely 15U, 
is completely reached in 0.8s. Finally, Figure 5 shows vc 


trend after a load change of 30% which occurs at 1.5s. 
As it can be seen, after an initial transient the system is 
stabilized again around the reference value. Even if the 
gains can be tuned to slightly change the performances 
achieved, in general it could be stated that the hyperbolic 
tangent function has better performances in the steady 
state, both as regards the disturbance response and the 
performance indicators (PF and THD). 

6 . PHASE INDEPENDENT CONTROL FOR 
MODULAR MULTILEVEL CONVERTER 

6.1 Model and Problem Formulation 

The MMC introduced by Prof. Marquart in Glinka and 
Marquardt (2003) has several advantages with respect to 
its predecessors, such as its high modularity, scalability 
and lower losses. It seems to be one of the most promising 
converters for bulk power transmission via HVDC links 
Ahmed et al. (2011). 

Nonetheless, controlling the MMC is no simple task. Sev¬ 
eral efforts have been oriented to propose suitable math¬ 
ematical dynamical models Antonopoulos et al. (2009); 



























Fig. 6. MMC Topology 

Harnefors et al. (2012), and control strategies for the 
converter under balanced operation. Some of the proposed 
control strategies have been implemented in dqo rotating 
reference frame such as Tu et al. (2011); Bergna et al. 
(2012). Nonetheless, it seems to be getting clear that there 
are significant disadvantages in applying dqo-hased control 
schemes for the control of the MMC differential currents 
since they do not facilitate complete phase-independent 
control of the converter state variables, hence the MMC 
potential will not be fully exploited Bergna et al. (2013). 
To overcome this drawback it is necessary to formulate the 
control problem in the abc frame—in which the MMC is 
controlled independently per phase—resulting in a track¬ 
ing problem instead of a regulation one. 


A = 


B, = 


be written in the form (1) with: 

0 0 0 

0 

0 0 0 0 

0 0 0 0 


■ 

rvda-\ 


2L 

,d = 

0 

0 


0 


0 fi; 


0 

-1 


^ 0 
,0 ^ 


0 417 


,Bo = 


0 

0 

1 

2C' 

0 


0 

-1 


- 1 ’ 

4L 


417 0 


Therefore, the matrix P satisfying Assumption (1) is 

r2L 0 0 O' 

0 L' 0 0 


P = 


0 0^0 


0 0 


0 ^ 
^ 2 . 


The passive output defined in (10) is 


y = 


xl Bj^ 


Px 


(29) 


2 yxii^XA - XiXA* + \x2*X3. - |a:2X3*_ 

The convergence condition (16) imposes that the rank of 
the matrix 


' C 
Q-- 


^^3* ^^4-*- ^^2-* 

2^"^* 4^3* 4^2* 2^1* 


Vm o_ 
0 Vi?' 
0 0 

0 0 


should be full. This is clearly the case if 

xl* - Tl V 0, 


(30) 


6.2 Generation of references 


The MMC converter studied in the paper is shown in Fig. 
6 . To obtain its mathematical model four state variables 
have been considered: 1) The differential or circulating 
current of the MMC idiff, 2) the grid or load current 3) 
the sum between the upper and lower capacitor voltages 
ucn ■= ucu + ucL, and 4) the difference between them 
MCA := Ucu — ucL- This yields the model 


R. 


’’-diff = 


Vdc 




R' 1 

jjiv - -^UCAU^ 


v; 


4L' 


UCY.UA 


pcc 

17~ 


1 1 

UCY =-^pdiffUY + -^PvUA 

UCA =^^'ivUY + -^idiffUA, (28) 


where the control signals are me and ma and the other 
signals and constants can be identified from Fig. 6 (Refer 
to Harnefors et al. (2012); Antonopoulos et al. (2009) for 
more details on the model). Also, henceforth Vpcc = 0 
since the system here studied is not considered to be 
connected to grid. In addition, notation for the equiva¬ 
lent capacitance C := C/N, the equivalent inductance 
L' := Ll2 Lioad and resistance R' := Rl2 + Rioad 
has been introduced for the sake of clarity. Defining x := 

[idiff iv ucY uca]^ and u := [me ma]^, the system can 


As indicated in Remark 1 obtaining the exact expressions 
for the desired trajectories cc* and their corresponding 
control signals m * that satisfy (2) is a very complicated 
task. Therefore, in this subsection several practical con¬ 
siderations have been made, which are widely adopted 
in the power electronics community, to approximate their 
solution. 

The reference equilibrium point for the grid current is 
imposed based on the desired voltage reference of the load. 
For the system under study with no grid connection (i.e, 
Vpcc = 0) such reference is found by means of simple 
phasor calculations; i.e., by dividing the internal e.m.f. 
voltage of the MMC (e„* = |e„*| by the equivalent 
impedance {Z' = R' + jwL'). 

X2* =-^J^^^=sin (wt - tan'T^^)) (31) 

where e,,* is the sinusoidal voltage that one wishes to 
apply to the load, with an amplitude that can vary from 
0 to The differential or circulating current reference 
Xi* is estimated by assuming low converter losses; hence, 
the mean value of the input ac power \Pac\ = |e«*X 2 *| is 
approximately equal to the mean value of the power at the 






































































dc terminals of the MMC \Pdc\ = VdcXi^. 
be expressed as 





Thus, a;i* can 


(32) 


Since xi* is constant in this case, the voltage that drives 
it may be expressed simply by Udiff* = -Rxi*. 


Also, it is possible to calculate the fluctuations of the sum 
ws* and difference wa* of the capacitive energy stored 
between the upper and lower arms of the MMC (see 
Harnefors et al. (2012) for more details on such equations). 
Thus, 

~ “f {Vdc ^"^dif ^1*) dt 

Jo 

Awa* = J {Vdc - 2udiff*) - 2e„*a;i*^ dt (33) 

Remark 5. A high-pass filter is needed for the term that 
is being integrated in equation (33) to leave out the 
power error caused by neglecting assuming no losses in 
the converter (32). 


The average value of rcs* and wa* is given as a reference 
from the user; typically, and 1 Ta 0 * = 0. 

Using the energy estimation, it is now possible to calculate 
the upper and lower arm voltages. 


N 

ucu* ~ \j [{Wy,4>* + Arcs*) + (1Ua0* + Awa*)] 

[n 

ucL* ~ \j + Awe*) — (1Ua0* + Awa*)] ( 34) 

The remaining state variables ue* and ma* may now be 
calculated as 


2^3* =UCU* + UCL* 


2 ^ 4 * =UCU* - UCL* (35) 

In addition, the upper and lower insertion indexes are 
calculated by 

2 ^v* Udiff* 

n«* =- 

'^CU'k 


ni^ = 


2 “1“ '^diffic 

Uclir 


(36) 


Finally, consistent with this methodology, the control at 
the desired trajectory is defined as: 


ni* — rill* F rii* 

U2* —Uuk* rii* (37) 


6.3 The Controller 


To summarize, the controller equations are: 


z = -y 

u = -KpU + KiZ + u* 


1 


y = 


2:i*2:3 2:4X3* F 2X2*X4 9 X 2 X 4 * 


with reference variables X 2 * defined in (31), xi* in (32) 
and X 3 *, X 4 * in (33), (34) and (35). Also, ui* and U 2 * are 
derived in (37) and (36). 


6.4 Simulations and experimental results 


The MMC simulation scenario has been set up in Mat- 
lab/Simulink using a high efficiency model Gnanarathna 


et al. (2011) to test the validity of the control. The con¬ 
siderations are the following: the converter has 2N = 10 
submodules, N in each arm (upper and lower). The input 
DC voltage is Vdc = 150U, the reference voltage e„* has 
an amplitude of ^ and the frequency set to 50Hz. The 
frequency of the balancing algorithm Antonopoulos et al. 
(2009); Harnefors et al. (2012); Glinka and Marquardt 
(2003) that balances the N capacitor voltages is set to 
20kHz. The internal capacitance, resistance and induc¬ 
tance are respectively set to C = 3.3mF, i? = 8D and 
L = lOmH. The load resistance and inductance values, 
respectively, are i? = 6D and L = 20mH. 

The figures presented show the controller performance for 
both the simulation and implementation scenarios. Figure 
7 shows the differential and load current, respectively, idiff 
and iy. The grid current tracks with significant accuracy its 
reference orbit. The differential current and the capacitor 
voltages achieve good performance; however, they are 
strongly influenced by high-order harmonic pollution. This 
is caused by the long switching dead-time requirement of 
the experimental setup. In Fig. 8 it is shown the values 
for the upper and lower capacitor voltages ucu and ucl- 
According to our calculations in the simulations, there 
is an error of around 3% between these voltages signals 
and their references. It can be concluded that is a direct 
consequence of the approximations made in the estimation 
process. Finally, in Fig. 9 are depicted the 2N voltages of 
the MMC capacitors while fig. 10 shows the output multi¬ 
level waveform of the converter. 

7. CONCLUSIONS 

In this article the trajectory tracking problem for power 
converters based on passivity foundations has been solved. 
Because of the use of passivity (of the incremental model) 
the controller is a simple PI. The stability results are 
global and hold for all positive definite gains of the PI. 
In fact, this outcome extends previous results obtained for 
the regulation case. 

The performance of the controller was tested by means of 
some realistic simulations and experiments for both the 
boost and the Modular Multilevel Converters models. A 
future work within this area of research includes the use 
of observers for partial state feedback and its succeeding 
application to power converters. 
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